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The liquid-gas phase transition and associated instability in two component systems are investi- 
gated using a mean field theory. The importance of the role of both the Coulomb force and symmetry 
energy terms are studied. The addition of the Coulomb terms bring asymmetry into a mean field 
and thus result in important differences with previous approaches which did not include such terms. 
, The Coulomb effects modify the chemical instability and mechanical instability domains shifting 

many features away from proton fraction point y = 1/2 to a value closer to the valley of (3 stability. 
These features are discussed in detail. 
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Thermodynamic properties of two and, more generally, multi-component systems are of interest in many areas of 
physics. For example, binary systems of two components are encountered in nuclei where the two components are 
neutrons and protons. While light nuclei have N = Z , heavy nuclei have N > Z and systems with large neutron 
excess such as neutron stars exist in astrophysics. In condensed matter physics two component systems appear in 
y—i ' liquid 3 He where the two components are spin-up and spin-down fluids. Such systems contain a more complex phase 
structure than one component systems. This letter is an extension of previous work concerning the nuclear case which 
was discussed from both the fragmentation point of view |i],|]j3) and from a mean field picture B . In the mean field 
picture we studied the phase surface associated with the liquid-gas phase transition incorporating the role of surface 
and Coulomb effects [Q. Previous mean studies H) °f two component systems were investigations of the role of the 
symmetry energy on the phase structure of the liquid-gas phase transition. While the symmetry energy plays an 
] important role in governing the proton/neutron ratios in the liquid and gas phases, favoring more symmetric systems 
• in the liquid phase than the gas phase, the effects of both the Coulomb interaction and surface energry were also 
shown to be important. This was first discussed in Ref. and more recently in Ref. ||. Here we will focus on the 
•t— j stability, metastability and instability of asymmetric systems of protons and neutrons. Such systems are of interest 
' for radioactive beams studies which can be used to explore nuclei at the limits of isospin asymmetry. In a recent 
article, V. Baran etal Q have pointed out some new effects associated with two component systems where a new 
kind of liquid-gas phase transition associated with chemical instability may appear. This new type of transition may 
manifest itself as an isospin distillation or fractionization as reported in Ref. [^],[)| where the gas phase is enriched 
in neutrons relative to the liquid phase. Such observations suggest that neutron diffusion occurs at the same rate 
or is more rapid than fragment production. These observations are based solely on symmetry energy effects, which 
make the liquid phase more symmetric, with the role of the Coulomb interaction not included. The importance of 
Coulomb effects driven by Z 2 terms versus symmetry energy effects arising from (TV — Z) 2 terms can be seen by just 
looking at the most stable nucleus for a given A. Using a liquid drop model, the Z with the smallest mass or greatest 
binding energy is given by Z = (A/2)/(l + B c A 2 ^/B s ) where B c is the Coulomb coefficient B C Z 2 /A^ 3 and B s is the 
symmetry energy coefficient B S (N — ZS 2 /A. When the Coulomb interaction is turned off Z is simply A/2. Moreover, 
the Coulomb interaction was shown pi to have a very dramatic effect on the specific heat associated with a first 
order liquid-gas phase transition, substantially reducing the peak previously reported in Ref. |3| where surface energy 
effects associated with one component systems played the dominant role. Moreover, the Coulomb force introduced 
another pair of coexistence densities having larger proton concentration in gas phase than in liquid phase due to its 
asymmetry in the mean field 0. 

In symmetric systems, mechanical instability from dP/dp < and chemical instability from dp /dp < at fixed T 
are simply related since 

Ndp = -SdT + VdP (1) 

where p is the chemical potential, S the entropy, P the pressure, N the number of particles and T the temperature. 
Thus pdp/dp = dP/dp at fixed T and the boundaries of both instabilities are common with dp/ dp = dP/dp = 0. 
The coexistence curve of the Maxwell pressure P versus density p or volume V is determined by the equality of the 
pressure and chemical potentials between the liquid and gas phase at the same T. In the region between the spinoidal 
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curve determined by the instability line and the coexistence line are metastable regions of supercooled liquid and 
superheated gas where nucleation processes can occur. These processes for one component systems were studied in 
Ref. , where nucleation processes gave a dynamical picture of the liquid-gas phase separation. Specifically classical 
nucleation theory |ll| looks at the change in the Gibbs free energy with and without a liquid drop made of A particles. 
This change is written as AG = (pi — p g )A + 4irR 2 a where pi is the liquid chemical potential of the drop, p g is the gas 
chemical potential from which the drop is formed, R the radius of the drop and a is the surface tension. A superheated 
gas has p g > pi , while equilibrium states have p g — pi . This change in Gibbs free energy has a maxima at a critical 
size drop whose radius is given by R c = 2a / (p g — pi)pi, where pi is the liquid drop density. The A = piV with V the 
volume of the drop. Drops with smaller R evaporate particles thereby lower their free energy, while drops with bigger 
R grow by accumulating particles and also lowering their free energy. Such processes in two component systems of 
neutrons and protons with additional Coulomb and symmetry energy effects have also been partly explored. Some 
preliminary studies of Coulomb effects in nucleation were done in Ref. @. For a two component drop, the AG is 



AG = (pi p - p gv )Z + (pi n - p gn )N + AttR 2 g + 5E C + 5E. 
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The SEc is the change in Coulomb energy with and without the drop and SE sym is the change in symmetry energy. 

In two component systems the coexistence line of a one component system now becomes a surface in (P, y, T), with 
y the proton fraction. This surface has several new elements associated with it that are not present in one component 
systems. For example, the isolated critical point of a one component system now becomes a line of critical points with 
y = 1/2 the critical point of a one component system. Associated with this surface is also a line of maximal asymmetry 
in the proton/neutron ratio and a line of equal concentration in liquid and gas phase The intersection of this 

(P, y, T) with a plane at fixed T leads to loops of P versus y as shown in Fig. 4 of Ref. pj . The chemical potentials 
p n and p p versus y at fixed T, P also have interesting behavior as shown in Fig. 3 of Ref. Q. For two component 
system, Eq.((l]) is replaced by 

y(dpp/dp) + (1 - y)(dp n /dp) = (\/p)(dP/dp) 
y(dp p /dy) + (1 - y)(dp n /dy) = (l/p)(dP/dy) 

In the first equation, y and T are held fixed, and in the second equation p and T are held fixed. We will use these 
figures to study the spinoidal instability region and its associated surface in (P,y 7 T). 

Some properties of the instability region can best be analyzed using some thermodynamic identities. We first note 
that the pressure and the chemical potentials are usually written as functions of p, y and T; P(p, y, T) and p q (p, y, T). 
One important identity arises from the equation of state P = P(p, y, T) which gives 



P,T 
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for the relation of y and p at fixed P and T. Since P, T, y and p are connected by an equation of state P(p, y, T), we 
can write 

dp q \ f dp q \ ( dp q \ f dp 



dy J PT \dy J p T ' \ dp J y T \dy J PT ^ 



Combining Eqs.(||) and (||), we can write (dP/dp) Vi T as 



dP \ _ - (dP/ dy)p yT (dpq/dp) y ,T 
d P J y.T {dp q /dy)p,T - (dp q /dy)p,T 



(4) 



Since (dp q /dy) p ,T is finite, singularites in (dp q /dy)p t T (chemical instability) give rise to zeroes in (dP/dp) y ,T- The 
zeroes of (dP / dp) V: T form the spinodal surface, boundary for mechanical instability region. In two component systems, 
mechanical and chemical instability regions are separated for asymmetric cases. Using Eq.(^), the chemical instability 
boundary ((dp q /dy)p t T = 0), determined by setting Eq.(|3|) equal to 0, is obtained from 

dP\ _ ( dP\ (dp q /dp)y t T ^ 



d P ) y ,T \dy ) p ,t {dp q /dy)p.T 

The following features should be noted. First, mechanical equilibrium ((dP/dp) y _T — 0) is determined by the condition 
that the left hand side of this equation is 0. Thus the mechanical and chemical instability boundaries are the same 
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when {dP/dy) p ^T = since (dfj, q /dy) Pt T ^ here. Due to Eq.(||), (dp/dy)p i r — too when (dP/dy) p .T = except 
on the mechanical instability boundary. This condition is also related to the equal concentration yE of liquid and 
gas phase on the coexistance surface [Q. For systems without Coulomb terms this occurs at y — 1/2. Secondly, for 
chemical equilibrium the right hand side can be obtained from the equation of state P — P(p, y, T) and the behavior 
of the proton or neutron chemical potentials p q — p q (p,y,T). The q will distinguish neutron results from proton 
results. We will use Eq.(||) to explore the boundary of chemical equilibrium. 

We begin our study by first considering a simpler system where Coulomb and surface effects are turned off. Using 
Eq.(54) of Ref. f| with C s — (surface term), C — (Coulomb term), x 3 = —1/2 and a = 1, then (dP/dp) y ,T is 
quadratic in p and is given by 



d PJy,T 
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The last term involving (2y — l) 2 does not contribute for symmetric systems. The terms involving A 3 arise from 
degeneracy corrections. The mechanical instability points can easily be obtaind by solving the quadratic equation 
(dP / dp) Vt T = to obtain the two densities that define the instability region for each T and y. These two roots become 
a single solution at a specific T c (y) at each y. This occurs when (d 2 P/dp 2 ) Vt T = 0. In one component systems and in 
symmetric systems with y — 1/2, this is the critical point T c . 




Density, p(fm 3 ) Density, p(fm 3 ) y 

FIG. 1. Chemical (thick solid curve) and mechanical (thin solid curve) instability curves at temperature T = 10 MeV with 
(upper boxes) and without (lower boxes) Coulomb effect of R = 8 fm uniform sphere. The dashed curves are the coexistance 
curve at T — 10 MeV. The dotted curves are for (djj, p /dp) Vt T = and the dash-dotted curves are for (dp n /dp)y,T = 0. The 
equal concentration yE = 0.4256 at T = 10 MeV. The curves are the constant T plane cuts of the corresponding surfaces in 
(p,y,T) space (left column), in (p,P,T) space (center column), and in (P,y,T) space (right column). 
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FIG. 2. The temperature dependence of the chemical and mechanical instability curves at fixed proton concentration y. The 
curves here are the constant y plane cuts of the corresponding surfaces in (p,y,T) space. The curves are the same as in Fig. 
[j] except the dashed curves. With the temperature dependent equal concentration Ve(T), all the curves coincide as shown by 
dashed line here (notice here that Ve(T) — 1/2 without Coulomb force). 

The calculations using Eqs.(53) and (54) of Ref. Q with degeneracy correction are shown in Fig. [I] and Fig. |^. 
The mean field force parameters are the same as used in Ref. Q except for using a larger radius R = 8 fm of a 
uniform sphere for a larger Coulomb force. These figures show that the mechanical instability boundary (thin solid 
line) is inside of the chemical instability region (thick solid line) and these coincide at the proton concentration of 
He = 0.4256 at T — 10 MeV (Fig.|l| and upper middle box of Fig.|J). The curves of (dp q /dp) y: T = also coincide with 
the instability boundary at ije- The coexistance boundary coincides with the chemical instability boundary at the 
line of critical points (peaks in Fig.|l| for a T — 10 MeV plane cut) but not with the mechanical instability boundary . 
We can see from Fig|2 that the chemical and mechanical instability boundaries coincide at the proton concentration 
Ue{T) which is on the line of equal concentration on the coexistance surface 0]. While all the curves coincide at for 
y = 1/2 (without Coulomb) as shown in the lower left box of Fig.|| the top middle box of Fig. with Coulomb exhibit 
a separation between these curves due to the temperature dependence of j/_e(T). All curves in the upper middle box 
of Fig. coincide at T = 10 MeV since Ue = 0.4256 at this temperature. From Figj^ and the left column of Fig.|l], we 
can see that all the four surfaces in (p,y,T) space coincide at the line of temperature dependent 2/e(T) and that the 
mechanical instability surface is inside of the chemical instability surface. Upper Fig. |^ and upper left box of Fig.|l| 
also shows that the role of (dp p /dp) y .T — and {dp n /dp) y .T — are exchanged when y crosses yE- 

To keep equations simple in the dicussions of various new features in two component system with Coulomb force, 
we will drop the degeneracy corrections. However, they are included in the actual calculations done and in the 
fiqures (Figs. |l| and ||) presented here. In this limit P = pT — ao(y)p 2 + 2a 3 p 3 where we have followed the notation 
of Ref. 1 10 1 with P now generalized to the case of asymmetric systems. The generalization changes do to ao(y) = 
a [l - (2/3)(l/2 + x ){2y - l) 2 ] with a = -(3/8)t and a 3 = i 3 /16. Then dP/dp = at 



P± 



(a (y) ± y/a? (y)-6a 3 Tj /(6o s ) (7) 



The two roots are equal at T c (y) — a^{y)/{Qa 3 ) and are p c {y) = P+ = p- = <Jo(y)/(6a3). At this point the pressure 
is P c (y) — ao(y)/(108 a i)- Under the same limits the chemical potentials are given by 

p q = Tln[(A 3 /7)(p/2 ± (2y - l)p/2)] - 2a oP ± (4/3) (1/2 + x )a a (2y - l)p + 3a 3 p 2 (8) 

with the upper sign for protons and the lower sign for neutrons and 7 = 2 for the spin up and down degeneracy. This 
equation for p q = p, q (p, y, T) when combined with the equation of state P(p, y, T) can be used to calculate the chemical 
potentials as a function p q = p q (P,y,T). Boundaries of the chemical instability region can also be determined using 
the condition dp q /dy = at constant T, P. Substituting the simplified equation of state and expression for the 
chemical potential into Eq.(|5|) leads to a cubic equation c 3 p 3 + C2P 2 + c\p + cq = with c 3 = 32(1/2 + xo)aoa 3 (l — y)y, 
c 2 = [6a 3 T — (32/3)(l/2 + £o) a o(l ~ v)v]-i c i = — (4/3)(l — x )aoT, Co = T 2 . Chemical instability occurs when three 
real roots exist, two of which will be postive, one negative. One negative root is present because co/c 3 is positive. 
At y = 1/2, the two positive solutions are the same as that for mechanical instability. When y is no longer 1/2, the 
mechanical instability and chemical instability regions are no longer the same. This result is shown in the lower part 
of Fig.|| for various choices of y with the degenercy correction of Ref. Q. It should also be noted that the coefficients 
of the cubic equation have a certain symmetry arising from the form in y which is y(l — y), i.e., interchanging y 
(protons) and 1 — y (neutrons) leaves the solution invariant (see curves for (dp q /dp) Vt T = in Figs, [l] and ^). This 
result is no longer true when Coulomb interactions are included. We now turn our attention to the simple EOS and 
p q , but now with Coulomb effects included in a simplified form which will illustrate its qualitative effects. 

When Coulomb effects are included in the approximation of a uniformly charged sphere of fixed radius R, the 
pressure equation contains an additional term Cy 2 p 2 and the proton chemical potential contains an extra term 2Cyp. 
The C, which depends on R, is taken to be a constant. For this choice, we can consider R as an effective range 
of the Coulomb force which is now independent of the density. This simplification reduces the complexity of many 
results. To make the Coulomb effects more visible in our figures (Figs. |and §), we use a large value of C = 231 
MeV/fm 3 which is appropriate for systems near U. These additions to P and p p change the equation for chemical and 
mechanical instability. Specifically the coefficients of the cubic equation for chemical instability for both protons and 
neutrons are now c 3 = [32(1/2 + x )a + 12C]a 3 y(l - y), c 2 = 6a 3 T- [(32/3)(l/2 + x )ag + (8/3)(l - x )a C]y(l - y), 
ci = [-(4/3)(l - x )a + 2Cy]T, c = T 2 . 

Some qualitative features of the effects of such a Coulomb behavior are the following (see Fig. |l| and Fig.|^ also). 
1) The mechanical and chemical instability curves are the same when (dP/dy) Pv T = which occurs at ys = 0.5/(1 + 
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C/[(8/3)(l/2 + Xo)ao]) and thus no longer at yE — 0.5 when (7^0. This equation for y has a similar structure to 
the equation that gives the Z which has the maximum binding energy at fixed A (valley of (3 stability) obtained from 
the liquid drop mass formulae or from Eq.(55) of Ref. 0). With higher order T dependence coming from degeneracy 
corrections, yE depends on T and yE — 0.4256 for xo — 1/2 at T = 10 McV. Explicitly, with first order degeneracy 
corrections, yE = 0.5/(1 + CJ [(8/3) (1/2 + x )a + TA 3 /(2v / 27)])- 2) For both larger and smaller values of y from yE, 
the chemical instability curve is outside the mechanical instability curve as shown in Figjl] and |^. The mechanical 
instability has its peak in T versus p at T c (y) = [do — (2/3)(l/2 + xo)aa(2y — l) 2 — Cy 2 } 2 / (6as) which is maximum 
at y = yE and coincide with chemical instability at yE- 3) The maximum critical T obtained from the peak in the 
chemical equilibrium surface of T versus p, occurs at y ^ 0.5, but is now shifted to y = ue(T). 4). At y = yE{T), 
the equilibrium curve or binodal curve of P versus y obtained by cutting the binodal surface at a given T meets at 
a point, i.e., the liquid and gas phase have the same proton concentration (see Fig. 4 of Ref. Q). For y > yE the 
gas phase is proton richer than the liquid phase when compared to a situation without Coulomb effects When 
Coulomb effects are turned off the binodal section meets at y = 1/2. Aty = 7AE^l/2 with Coulomb effects included 
a two component system will behave as a one component system, keeping the proton/neutron ratio the same in both 
liquid and gas phases. 

In summary, we have presented calculations based on a mean field theory with asymmetry which include the role of 
Coulomb force on the liquid-gas phase instability regions for both the mechanical (spinodal) surface and the chemical 
(diffusive) surface. The Coulomb force modifies these surfaces so that their peaks are no longer at proton fraction 
y = 1/2 but shifted to a smaller value called yE which is close to the proton fraction in the valley of /3-stability. The 
mechanical instability surface is always inside the chemical instability surface and the two surfaces are the same at 
the equal concentration yE (allowing for T dependence in ys). At yE, the peak in the chemical instability surface 
and binodal equilibrium surface also coincide. For systems with y ^ yE, the peak in the T versus p curve for chemical 
instability and chemical equilibrium - binodal or coexistance curve - coincide, but the peak for mechanical or spinodal 
instability is below that of the other two. At y = yE, a phase transition in a two component system behaves the same 
as a one-component system, keeping the proton fraction at yE in both phases during the entire transformation. 
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